Anomalous thermodynamics of Coulomb interacting massless Dirac fermions in two 

spatial dimensions 
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It is argued that the specific heat of N massless Dirac fermions in 2 spatial dimensions inter- 
acting with 1/r Coulomb interactions is suppressed logarithmically relative to its non-interacting 
counterpart. The (dimensionless) coefficient of the logarithm is calculated in a closed form in the 
leading order in large N expansion, but to all orders in the effective fine structure constant, hf, a 
procedure which takes into account finite temperature screening. This effect is expected to occur in 
a single layer graphene embedded in a dielectric medium. Its dependence on the dielectric constant 
is calculated analytically. 
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Understanding properties of quantum matter confined 
to two spatial dimensions has been at the forefront of 
theoretical physics [l| . The lowered dimensionality is be- 
lieved to increase the role of interactions and, together 
with quantum statistics, lead to anomalies in various 
physical observables 0, d, 0| • 

In the case of Fermi-Dirac systems, the role of (short 
range) interactions has been studied in the pioneering 
work of Landau[5] . What is now understood as the Fermi 
liquid paradigm explains the low temperature proper- 
ties of the degenerate Fermi-Dirac systems as arising 
from a dilute gas of weakly interacting quasiparticles 
whose quantum numbers are the same as those of non- 
interacting fermions. In this context, it is also well known 
that a degenerate assembly of electrons with a rotation- 
ally invariant Fermi surface, in either three or two spa- 
tial dimensions, interacting with 1/r Coulomb interac- 
tions, can be described by the Fermi liquid theory Q. 
Ultimately, this is due to the ability of the mobile carri- 
ers to screen the long ranged interaction beyond the dis- 
tances of the order of the inter-electron separation. One 
of the most direct observable consequences for a metal is 
the linear temperature dependence of the specific heat at 
asymptotically low temperatures i.e. limj^o cv/T = 7 
0,|. 

The situation is less clear in the case of a semimetal, 
which can be thought of loosely as the Fermi surface 
shrunk to a point. More precisely, in this case, a con- 
duction band and a valence band typically touch at a 
discrete set of points with dispersion which vanishes lin- 
early with the wavenumber near each such point, being in 
a sense a critical point between a metal and an insulator. 
If the chemical potential lies at such a point of degener- 
acy, the absence of the Fermi surface implies that, unlike 
in a metal, the long ranged interactions are screened only 
by thermally excited carriers, with the screening length 
of the order of the thermal length, Hvf /ksT. For a typ- 
ical Fermi velocity ~ 10 6 m/s, even at room temperature 
the screening length can therefore be much longer than 
the interatomic spacing. The effect of the Coulomb in- 



teractions on 3D semimetals was investigated in Ref. [9( 
where it was argued that, unlike in a metal, the Coulomb 
interaction causes the energy spectrum to differ from a 
purely linear spectrum by a logarithmic factor. 

In this work, I study the effects of the Coulomb in- 
teractions on the thermodynamic properties of massless 
Dirac fermions in 2D, which can be used as a low en- 
ergy description of a single-layer graphene, and find that 
the poorly screened electron-electron interaction leads to 
anomalies in the quasiparticle thermodynamics. 

In the absence of interactions, the free energy density 
is (up to a temperature independent constant) 
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where N is the number of the 2 component Dirac fla- 
vors (N = 4 in the single layer graphene) and C( n ) is t ne 
Riemann zeta function. The strength of the electron- 
electron interaction is given by the dimensionless cou- 
pling constant of = e 2 /(zTivf), which is the effective 
fine structure constant, where e describes the polarizabil- 
ity of the surrounding medium. If e is large, then olf is 
small, and to first order in ap , vf receives a logarithmic 
correction jVpotF m []nj where A is a short wavelength 
cuttoff. Such logarithmic enhancement of the otherwise 
linear dispersion near the Dirac point, would imply log- 
arithmic suppression of the specific heat. It is not clear, 
however, whether such simpleminded reasoning suffices. 
Similar reasoning in the case of an electron gas with a 
Fermi surface would also lead to a logarithmic suppres- 
sion of the effective mass, the result known to be incorrect 
due to screening. Therefore, the key issue in this context 
is the effect of the thermal screening on the specific heat 
in the 2D Dirac case. 

To answer this question, I set up a large N expansion 
while keeping cyfN fixed and arbitrary. To the leading 
order in large N, this is equivalent to a random phase 
approximation (RPA). I find that the free energy density 



receives a logarithmic correction 
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where A = apNTr/8, the large temperature cutoff is re- 
lated to the short wavelength cutoff by HvfA — ksTuy , 
and the function g is given by the Eq. (fl~2")) . 

To the same order in the large N expansion, the above 
result can be interpreted as a free energy density of 2N 
modes with dispersion 
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I find that as A — > 0, 77 = 2X/(wN) while for A — > 00, 
7/ = 8/(tt 2 N) (see Fig[T]). The expression (JTJ can be 
understood as the first non-trivial term in the Taylor ex- 
pansion of 
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Note that in the limit of A — > 0, the above expres- 
sion coincides with the free energy calculated within the 
Hartree-Fock approximation. As T — > 0, the free energy 
in (J2J vanishes as -T 3 / In 2 t[Ii|, E3- Thus the suppres- 
sion of the specific heat cy = —Td 2 f/dT 2 relative to 
the non-interacting case persists when the polarization 
effects are included. 

The logarithmic divergence of the particle's group ve- 
locity at low wavenumbers may seem to violate causality. 
This is an artifact of the approximation which treats the 
1/r 2 Coulomb forces as instantaneous. Were we to in- 
clude the retardation effects introduced by coupling the 
fermions to a (quantum) 3D electromagnetic gauge field, 
the Fermi velocity would not grow without bound at low 
k, but instead saturate to the speed of light [HI]. How- 
ever, since the growth of Dp is only logarithmic, the ef- 
fects of retardation on the specific heat would be practi- 
cally unobservable for c> vf- For all practical purposes, 
the above non-relativistic expression will therefore span 
the physically relevant regime. 

Below, I justify the above claims. For convenience, 
I work in units where h — k B = vp = 1 and restore 
the physical units in the final expressions. Then, the 
Hamiltonian is 



H = Ho + V. 
The free part of H is 
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FIG. 1: The coefficient 77 = i [£\ 2 g(\) - ±A of the log- 
arithmic increase of the free energy in Eq. ^ vs. A = 

„2 , 



iVe*/16 = Nj € jr v . The closed form expression for g is given 
in Eq.jnj). For small A, 1) = 2X/(nN) while for A -> 00, 
77 = 8/(ty 2 N). 



where p = — iV, er's are the Pauli matrices, and ipj(r)'s 
are two component anticommuting Fermi fields. N is 
the number of the Dirac fermion flavors which, in the 
condensed matter setting, equals the number of valleys 
times the number of spin directions (N — 4 for a single- 
layer graphene) . The interaction part of H comes from 
the (3D) Coulomb interaction between the charge density 
fluctuations and reads 
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Here Sn(r) = J2j' , Pj{ r )' t Pj{ r )i an d e is dielectric constant 
of any potential surrounding insulator. Note that in our 
units e 2 /e, which is the analog of a fine structure con- 
stant, is dimensionless. 

Under a scale transformation r — > br the Hamiltonian 
(J3j> appears to transform as H — > 6~ 1 7i. This might sug- 
gest that even in the interacting theory, the temperature 
dependent part of the free energy density / goes as T 3 . 
However, since the theory is well defined in the UV only 
in the presence a short distance cutoff A -1 , this symme- 
try and the T 3 form will in general be violated. 

In the large N expansion the leading non-trivial cor- 
rection to the free energy density due to interactions is 
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where e 2 = 27re 2 /e. njj e< (g, uj,T) is the retarded polar- 
ization function of the free fermions. 



For n 2 > q 2 : 
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while for q 2 > il 2 
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At t = 0, n r Q et ( q ,n,o) 
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Note that the interaction correction to the free energy 
has the scaling form 

Sf(T, N, e 2 , A) = T 3 T(e 2 N, A/T), 

and in the low temperature limit of interest ^ ^> 1. I 
now proceed with the evaluation of Eq.© in this limit. 

Consider first the static regime q 2 > SI 2 . As can be 
readily seen from Eq. (fl"0|) , for large q QmU^ et (q,i},T) is 
small compared to q/e 2 and the real part, and the latter 
can be approximated by 3?eilQ e *((7, ft, 0) to the same or- 
der. As a result, the contribution to Sf from the static 
regime is 
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where the last line represents only the most divergent 
contribution. 

Similarly, the most singular contribution from the dy- 
namic regime, Q, 2 > q 2 , can be evaluated by expand- 
ing (J6j) to first order in deviation of 3mIiQ e *(q, fi, T) 
from 3mriQ e *(q, ft, 0) as well as to the first order in 
5ienQ e *(<7, Q, T), both of which vanish at large q. The 
algebra is somewhat tedious and I just state the final 



result for the asymptotic expansion in A/T: 
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where the function g(x) is 
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This integral can be evaluated by a substitution followed 
by a partial fraction decomposition and one finds 
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Combining everything together I find that up to O(-k) 
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where A = Ae 2 /16. This is the result quoted in Eq.(TT]). 
The corresponding function jj = i [i|A 2 .g(A) — -4-1 is 
plotted in Fig. (fT]). Note that the contribution from 
small q, which also includes the effect of the thermo- 
plasma mode [13], gives only terms of order T 3 . 

If we interpret the above expression as an expansion 
of the free energy given in Eq. ((2|) , the low temperature 
specific heat is subquadratic. Fig. |2]) shows this sup- 
pression for realistic material parameters pertinent to the 
single-layer graphene. 
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ability of the charge carriers to screen the long range 
interactions, and as a result the density of states is effec- 
tively suppressed near the Dirac point. The experimental 
observation of the effect would be an important step to- 
wards our understanding of critical massless matter. 

I wish to thank Professor Tesanovic for useful discus- 
sions. 



FIG. 2: The suppression of the specific heat (per Dirac cone) 
relative to the non-interacting result Cy = 9 1^ k%T 2 / (Hvf) 2 
(dashed line) vs. T/T uv . This plot is for^JV = 4, v F = 
10 6 m/s, e = 1 which gives r\ ~ 0.14. 



It is interesting to make a connection between this re- 
sult and a seemingly unrelated one of Ref.#[l3|, where 
the authors studied the specific heat of N massless 2D 
Dirac fermions coupled only to the space component of a 
2D gauge field within the large N expansion. Unlike our 
Eq.([T]), the correction to the free energy found numeri- 
cally in [ill has the opposite sign, i.e. the specific heat 
is enhanced relative to the non-interacting case. The 
numerical coefficient of the T 2 In ^ enhancement of cy 
is 2.79 15], while in this work I find that upon setting 
A — > oo the coefficient of the specific heat T 2 In T suppres- 
sion is 72£(3)/7r 3 f» 2.791. The opposite signs of these 
terms can be understood as follows: setting A to infinity 
in Eq.#|JT]), recovers the leading logarithmic contribution 
to the specific heat of QED3 within the large N expan- 
sion arising from coupling solely to the time component 
of the gauge field. However, due to the Lorenz invariance 
of QED 3 , the leading logarithmic contributions from the 
time and the space components of the gauge field must 
cancel exactly 16j . It is in this sense that the result in 
Eq.(p} is consistent with Ref.flq]. 

In conclusion, I found that the thermodynamic signa- 
ture of the long range electron-electron interactions in a 
two dimensional semimetal are much more pronounced 
than in the case of a metal. This stems from the in- 
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